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Abstract 



We study the Euclidean two-point correlation function G q (x) of the topologi- 
cal charge density in QCD. A general statement based on reflection positivity 
tells us that G q {x) < for x / 0. On the other hand the topological sus- 
ceptibility \q = / d d xG g (x) is a positive quantity. This indicates that G q (x) 
developes a positive contact term at x = 0, that contributes to the determi- 
nation of the physical value of Xq- We show explicitly these features of G q (x) 
in a solvable nontrivial continuum model, the two-dimensional CP N ~ l model 
in the large- iV limit. A similar analysis is done on the lattice. 
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I. INTRODUCTION 



By the axial anomaly equation, matrix elements and correlation functions involving the 
topological charge density operator can be related to relevant quantities of the hadronic 
phenomenology. The most important phenomenon related to the topological properties of 
QCD is the explicit breaking of the U(l) axial symmetry, with the result that the theory 
contains neither a conserved U(l) axial quantum number, nor an extra Goldstone boson. 

In Euclidean space the topological charge density is 

g 2 

Q&) = T^e^Tr F„ u {x)F p(J {x). (1) 

The correlation function 

G q (x - y) = (q(x)q(y)) (2) 

is important for understanding the breaking of the U(l) axial symmetry. According to 
Witten's argument |l| the breaking of the U(l) axial symmetry should occur at the lowest 
non-planar level, i.e. at the next-to-leading order of the 1/N expansion, and require a 
nonzero large- N limit of the topological susceptibility Q 

X g = J d A x G q (x) (3) 

in the pure SU(N) gauge theory (see also the recent Ref. ||). An extension of this idea 
relates Xq °f the large- N pure gauge SU(N) theory to the rj' mass |4]]. Many numerical 
studies, based on Monte Carlo simulations of lattice formulations of the pure gauge theory, 
have been devoted to the estimate of Xq ( see e -g- Refs. |]5|-[T6|). 

Arguments based on reflection positivity [T| tell us that G q (x) < for \x\ > (this 



was already noted in Refs. jTSjOJ). On the other hand Xq > trivially from its definition. 
These facts indicate that there is a positive contact term at x = 0, that contributes to 
determine the physical quantity Xq- 111 order to investigate this issue, we consider the two- 



dimensional CP N 1 model p0,21|, which is an interesting theoretical laboratory for studying 



general topological properties. We will calculate G q (x) in the large- N limit, providing an 
explicit analytical example where the main properties that should characterize G q (x) in 
QCD are realized. As argued for QCD, G q (x) presents a short- distance singular behavior 
characterized by a positive contact term at x = and a negative diverging approach for 
\x\ —* 0. Neverthless, the low-momentum behavior of G q (x) and in particular its moments 



Xq ,j = j d d x (x 2 yG q (x) (4) 
(where Xq,o = Xq) turn out to be well defined and finite. 



1 The physically relevant topological susceptibility is correctly defined by Eq. (fj|) in terms of the 
Euclidean correlation function G q {x) (see e.g. Ref. Q and references therein). 
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The paper is organized as follows. In Sec. we discuss some features of G q (x) in the 
QCD theory, that can be inferred using reflection positivity, and other general arguments 
such as perturbation theory and renormalization group. In Sec. [TTT]we compute the large- N 
limit of (q(x)q(y)) in the two-dimensional CP N ~ X models. Its main features are discussed 
and compared with those of the correlation function G q (x) of QCD. In Sec. [IV] we study the 
continuum limit of correlation functions of lattice discretizations of the topological charge 
density operator. An explicit calculation is presented in the large- N limit of a lattice for- 
mulation of two-dimensional CP N ~ X models. 



II. REFLECTION POSITIVITY 

In order that Euclidean correlation functions can be continued back to Minkowski space, 



they have to obey a positivity condition: the so-called reflection positivity |fL7| , p2|| . The 
general statement concerning reflection positivity is that 

((QF)F) > 0, (5) 

where is the antilinear reflection operator consisting in an Euclidean time reflection and 
a complex conjugation, and F is an arbitrary gauge invariant function of the fields having 
support only at positive Euclidean times (see also Ref. f23j). As a consequence of the intrinsic 
odd parity of q(x) under reflection, 

Qq(x) = -q(8x) = -q(x 1 ,x 2 , x 3 , -x 4 ), (6) 

reflection positivity states that 

G q {x) < for \x\ > 0. (7) 

This fact holds for any operator that is intrisically odd with respect to reflection symmetry 
in the Euclidean space. 

The asymptotic large- and short- distance behaviors of G q (x) can be inferred by general 
arguments. At large distance G q (x) should decay exponentially (in the presence of fermions): 
G q (x) ~ e~"V r apart from negative powers of r = \x\. Simple dimensional, perturbative 
and renormalization group arguments tell us that for r —>■ 



G q (x) . . 

r 8 (mr) 



i o[-L 

mr 



(8) 



where c is a negative constant. The logarithms can be related to a running coupling constant, 
indeed in perturbation theory G q (x) is 0(g 2 ). Since the topological susceptibility Xq is 
positive (x q = in the presence of a massless fermion) and G q (x) < for x ^ 0, G q (x) should 
develope a positive contact term at x = 0, that compensates the negative contribution of its 
integral for i ^ and makes Xq positive. Inspite of this singular short- distance behavior, 
the low-momentum behavior and in particular the moments Xq,j °f G q (x), are conjectured 
to be well defined and finite. In the next section we will see in a solvable model how these 
features can coexist. 
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In the lattice formulation of the theory one may define two versions of the reflection 
symmetry that are equivalent in the continuum limit: site- and link-reflection simmetry. 
We recall that reflection positivity is essential on the lattice for the existence of a self- 
adjoint Hamiltonian at finite lattice spacing defined from the transfer matrix. Site- and 
link-reflection symmetries are both satisfied by the Wilson lattice action of SU(N) gauge 



theories |22 ], and by Wilson fermions with Wilson parameter r = 1 \24\. So Eq. (|5|) must 
hold also on the lattice for finite lattice spacing. A lattice discretization of the topological 
charge density operator is for example ||: 

1 ±4 

Ql{x) = - 4 2 e^ pa TT[U^U P a], (9) 

pvpa=±l 

where 11^ (x) = + jj^U^x + v)Ul(x) is the plaquette operator and U^(x) is the 

link variable, the sum is done over positive and negative directions. One may easily verify 
that qi{x) trasforms as q(x), cf. Eq.@, under site-reflection. Again reflection positivity 
tells us that G q (x) = {qL{x)qL{fy) is negative for x ^ (at least along the directions of the 
lattice and when there is no overlap among the link variables of the two operators). In the 
continuum limit the lattice correlation G^(x) should reproduce the continuum correlation 
function G q (z). This will be discussed later, in Sec. fV[ 



III. THE LARGE- AT LIMIT OF {q(x)q(y)) IN TWO-DIMENSIONAL CP^ 1 

MODELS. 

The arguments of the previous section can be also applied to the two-dimensional CP 1 ^^ 1 
models as well. In the following we will consider their large- A limit and will calculate the 
leading order of the two-point correlation function of the corresponding topological charge 
density (i.e. 0(1/A) ). We will show explicitly that G q (x) develops a singular behavior 
at the origin consistently with the reflection positivity requirement G q (x) < for x ^ 0, 
and the positivity of the corresponding topological susceptibility, i.e. of its space integral. 
Nevertheless, the low-momentum behavior of the G q (x) turns out to be well defined without 
the need of special subtractions. This provides an explicit example where the conjectured 
main features of the two-point function G q (x) of QCD are verified. 

(jpN-i moc iels are defined by the action [pO,pT 



N 



S = —J d'xD^zD^z, (10) 

where z is a A-component complex scalar field subject to the constraint zz = 1, and the 
covariant derivative = + iA^ is defined in terms of the composite gauge field = 
izd^z. Such models present a U(l) gauge invariance related to the local transformations: 
z(x) — > e ia ^z(x) and A^(x) — > A^(x) — d^a(x). 

Two-dimensional CP N ~ l models are interesting because they present many features that 
hold in QCD: asymptotic freedom, gauge invariance, existence of a confining potential be- 
tween non gauge invariant states (that should get eventually screened by the dynamical 
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constituents), and non-trivial topological structure (instantons, anomalies, 6 vacua). More- 
over, a pleasant feature of these models is the possibility of performing a systematic 1/N 
expansion around the large- N saddle point solution pO| , |2T|j25|] . 

In two-dimensional CP N ~ X models one defines the topological charge density operator 

i i 

q(x) = — e Mi/ D^zD u z = — d^A u , (11) 

and its two-point correlation function G q (x—y) = (q{x)q(y)) . q(x) is a renormalization group 
invariant operator, i.e. it does not have an anomalous dimension. Like the topological charge 
density of QCD, q(x) transforms as 

Qq(x) = -q(9x) = —q(xi, —x 2 ) (12) 

under reflection symmetry. Therefore, as a consequence of reflection positivity, G q (x) < 
for i^O. 

Let us introduce the second-moment correlation length £ associated with the two-point 
correlation function of the operator Pij(x) = Zi(x)zj(x), 

G P {x-y) = (TrP(x)P{0)), (13) 



Jd 2 x \x 2 G P {x) 
J d?x Gp(x) 



(14) 



We will use £ as length scale in the following calculations 0. 

The Fourier transform of G q (x) can be written in terms of the propagator of the composite 
field A^P^ip), 



(15) 



Substituting the leading large- iV expression of P* P,PJ25| , one finds 



NG q (p) = —p 



2n 



u(p) In 



u{p) + 1 
u{p) — 1 



u{p) 



3pX 



2 £2- 



(16) 



It is worth noting that in two-dimensional CP N 1 models the two-point correlation function 
of q(x) is also related to the correlation of two Wilson loops constructed with the abelian 



field in the limit of small area p5 



2 The second-moment correlation length of Gp(x) turns out to be more suitable for a 1/N- 
expansion than the mass-scale determined from the large-distance exponential decay of Gp(x), 
due to its analytical properties in 1 /N ]25| . 
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Since it is convenient to work with dimensionless quantities, we define 

B(k)=Z 2 NG q (jp = k/Z), 
i.e. we use £ as unit of length. B(k) has the following asymptotic behaviors 



B(k) 



k 2 



B(k) = - 



27rln(6P/e 2 ; 
1 



+ 



( 1 



\ln k 



2tt IOtt 



27 
350tt' 



0(k e ), 



(17) 

(18) 
(19) 



for large and small momentum respectively. 

The singular behavior of G q (x) at small distance is already apparent from the asymptotic 
behavior (pj|) of its Fourier transform. The calculation of the large- N limit of G q (x) requires 
to perform the Fourier transform of the expression (|T6|). As before let us introduce the 
dimensionless quantity: 



c(x) = ecM) 



d 2 k 



Ak-x 



B(k). 



The moments of C(x) 



X.. 



g<i 



d 2 x (x 2 yc(x) = e (1 - j) x q>j 



(20) 



(21) 



(where Xgj = / d 2 x (x 2 yG q (x)), and therefore of G q (x), can be easily obtained from the 
expansion of B(k) in powers of k 2 , cf. Eq. (|TT?D. In the large- N limit one finds 



Xq,l 



1 



2nN 
6 

~5^N 



1 

N 2 



O 



+ oi W 2 



(22) 
(23) 



etc... We mention that Xg an d Xg,i are known to 0(1/N 2 ) P^J2^7[] . 



By rotational invariance C{x) depends only on r 



\x\: 



C(x) 



1 

27 



kdkJ (kr)B(k). 



(24) 



In order to evaluate C(x) for finite r ee |x|, it is convenient to modify the integration contour 
in the complex k plane, moving it along the imaginary axis (when r > 0) [^J. For r > 
one can then write 



C(x) 



2tt 2 



dtK (tr)t 3 v(t) 



v(t) In 



l+v(t) 
1 - v(t) 



+ TX 2 v{tf 



(25) 



where v(t) — y 1 — Since -ft'o( ::c ) > 0, Eq. ( pop shows that C(a;) < for r > as expected. 
In Fig. HI we show the function C(x). The integral representation ( p5|) for C(x) holds only 
for r > 0. For r = the contour rotation leading to the integral representation 



misses 
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the contribution from the path at infinite distance, that is not suppressed anymore and that 
should represent the positive contact term. 

The integral representation ( p5|) allows us to derive the asymptotic behaviors of C(x). 
At large distance C (x) decays exponentially 0: 



C{x) 



1 e" 



24tt r 2 



l + 0(- 

r 



For r — > 0, C(x) diverges as 

C(x) 



2tt 2 r 4 (In-, 



1 + 



lnr 



(27) 



(28) 



One can infer this short-distance behavior also by calculating the leading order of perturba- 
tion theory that is given by 



2vr 4 r 4 



(29) 



Then, using renormalization group arguments, one replaces the coupling g with a running 
coupling constant g(r) « 7r/ln(A/r). 

The diverging negative integral of C(x) for r > must be compensated by a diverging 
positive contribution of the contact term at r = (that got lost in the integral contour 
rotation performed to evaluate C(r) for r > 0) so that 



d 2 x C(x) 



(30) 



The integrals for r > r > are finite; for example numerically one finds J r>1 d 2 x C(x) = 
—0.0503..., f r>1 / 2 d 2 x C(x) = —0.1653..., etc... Notice that as a consequence of the r~ 4 
short- distance behavior and the positivity of Xqi we have formally f\ x \ < $d 2 x C(x) — >■ +oo 
for 5 — > 0. Thus a 5- like distribution cannot represent the contact term. The behavior at 
x = should be described by more complicated distributions acting in a finite interval Q. 
Similar considerations hold in QCD. 



3 We used the asymptotic behavior 



TrV/ 2 



- 



2xJ 



K (x) : 
for x — > oo. 
4 An example of such distributions may be 



l + O 



(26) 



lim 



P £ (d)S(x)-X eT1I f (In |x|) 0(x-e) 



where the polynomial P £ {d) and A e are appropriate functions of e and the limit e 
considered in a weak sense, i.e. after performing the integral with the test function. 



(31) 
must be 
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For comparison one may calculate the large- N limit of Gp(x), cf. Eq.fllB]). Since the op- 
erator Pij(x) is not renormalization group invariant, in order to construct a renormalization 
group invariant function, one may define 

/W-^ (32, 

where xp — Gp(0) = ^g 2 (, 2 is the corresponding susceptibility, and therefore XpH 2 is the 
renormalization constant at zero momentum of the operator Pij(x). In the large- iV limit 
one finds (in units of £) 



whose asymptotic behaviors are 



ln(6/c 



™ = -$r +0 {-w)- (34) 

f(k) = 1 - k 2 + 0(k 4 ) (35) 
for large and small momentum respectively. After some manipulations, for r = \x\ > 

m = -[Z (36) 

JK ' tv(t) v ; 

(where again an integral contour rotation of the Fourier transform of f(k) has been per- 
formed), with the following asymptotic behaviors 

l„-i/fr 



f(x) ~ VVf for r -f oo, (37) 

/(x)~(lnr) 2 for r -> 0. (38) 

In the integral contour rotation the contribution to the contact term of the path at infinite 
distance is again lost. Comparing with the behavior of the topological charge density cor- 
relation function, we note that the singularity at r = is much softer in this case. This 
is already apparent by comparing the large momentum behaviors of f(k) and B(k), cf. 
Eqs. fl3|) and ©. 

IV. THE TOPOLOGICAL CHARGE DENSITY CORRELATION FUNCTION ON 

THE LATTICE. 

In a pure gauge theory q{x) is a renormalization group invariant operator, thus G q (x) 
is the universal continuum limit of corresponding lattice correlations (qL(x)qz,(y)) , when 
expressed in the appropriate units. In the presence of fermions a nontrivial anomalous 
dimension is generated due to the mixing with the operator d^j^x) ||28|| . One may still 



define a universal RG invariant function from the topological charge density correlation 
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function, using the fact that the anomaly equation (in the chiral limit) d^j'. 



i2N f q{ 



x 



does not get renormalized 

Further subtleties are present on the lattice. For example let us consider the case of a 
pure gauge theory. A more careful analysis of the continuum limit of the matrix elements 
of a generic lattice discretization qh{x) of q(x), such as Eq. (|9|), leads to the relation UTOf 

q L (x) = a*Z(g)q(x) + 0(a 6 ), 

where a is the lattice spacing, g is the bare lattice coupling and Z(g) 
renormalization function. Thus we have 



(39) 

1 + 0(g 2 ) a finite 



G L Ax) = a 8 Z(g) 2 G q {ax) 



l + 0(a 2 



for ax 7^ 0. 



(40) 



In order to evaluate the moments of G q (x) and in particular the topological susceptibility, 
one also needs the contribution at x — 0. This is strongly affected by lattice artifacts that 
eventually become dominant in the continuum limit. So a careful subtraction is required 
(see e.g. Ref. ||3T)| , |il"|j ). In order to overcome these problems, geometrical definitions have 
been proposed |31|j6|], that represent appropriate interpolations among the discrete lattice 
variables, and are in general non-analytical and non-single-valued functions of the lattice 
variables. These lattice estimators are not affected by renormalizations. On the other hand 
the corresponding topological susceptibility, and probably also higher moments, turns out 
to be sensitive to unphysical lattice defects [52], that may spoil their continuum limit in 
some cases (see also Refs. fl33| 



In order to investigate how G q (x) is recovered in the continuum limit of lattice corre- 
lations Gq(x — y) = (qL(x)qL{y)) > we consider a lattice formulation of the two-dimensional 
(jpN-i moc [ e is. We consider the lattice action |35L 



S, 



(41) 



where /3 = 1/g, z n is a iV-component complex vector, constrained by the condition z n z n = 1, 
and A„ M is a U(l) variable. One can easily prove that site- and link-reflection positivity holds 
for the lattice action ([H]). The lattice formulation ( |4~1~D of CP N ~ X models turns out to be 
particularly convenient for a large- N expansion [0^ - ^O| , |35|| . 

In an infinite lattice (free boundary conditions are assumed) one may consider the fol- 
lowing discretization of the topological charge density operator 



47T 



(42) 



where 6 n ^ is the phase of the field A„ )M , i.e. X n ^ = e Using the property of 9x(n) under 
site-reflection one can prove that G^(x) = (qz(x)qL(0)) must be negative for \x\ > (at 
least along the directions of the lattice). At large iV one can explicitly show that (^(n) has 
the correct continuum limit, and no lattice renormalizations are necessary [3S[35|. Thus in 
the continuum limit and for \x\ > one expects 



C L {x) 



?NG L AxZ) 



c(*) + o(r 2 ), 



(43) 
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where C(x) is the continuum function defined in Eq. (|20|), and £ the second-moment corre- 
lation length associated with the lattice correlation function Gp(x) = (Tr P(x)P(0)). 
Gq(x) can be calculated in the large- N limit. Straightforward calculations lead to 



NG L q (k) =Y,e ikx NGl 



x 



(27T) 



(44) 



where |37|-|39]j 



* d 2 q 2J2^cosq fl f« d 2 q 4 ^ sin 2 + fc M /2) 
-7r (27r) 2 g 2 + m 2 . 



(2tt) 



(q + k) + ml 



and q 2 = J2u Qu = sin 2 ^^). The parameter m 2 , is related to (3 by 



(45) 



P 



* d 2 q 



K 



(2n) 2 q 2 + m 2 2vr(l + mg/4) \l + m 2 /4: 



(46) 



(K is the standard elliptic function), and to the second- moment correlation length associated 
with Gp(x) by 



2G£(0) dk 



92 G L P (k)\ k =o = -A(a)(0)^A^(fc)| 



fc=0; 



(47) 



/7T 
-7T 



d 2 g 



1 



(2vr) 2 [ ^ + m 2 ] 



(q + k) + ml 



G%(x) can be evaluated by Fourier transforming G q (k). The lattice topological susceptibility 
is obtained by summing G q (x) over the lattice sites, 



X L q =Y,G L q (x) = G L q (U). 

X 



(48) 



An analysis of the above large- N expressions leads to the following main results. 

(i) The continuum limit of £, 2 G L (k) at A;£ fixed is B(k£), cf. Eq. flT7|). Indeed at large £ 



eG L q (k)=B{ki) + 0{C 2 ). 



(49) 



This may be seen by performing an asymptotic expansion of G q (k) (at fixed k£) in powers 
of £~ 2 , following the procedure outlined in Ref. f4"0]] . 

(ii) Gg(x) is negative everywhere for x ^ 0, consistently with reflection positivity. 

(iii) At fixed physical distance r = x/£ > the continuum limit exists and it is given by 
C(r), in agreement with Eq. flUf ). Notice that the convergence is not uniform in r. Moreover, 
from Eq. ( fj9|) it follows that 



(50) 



where 3 - are the lattice moments of G q (x) and x g , 7 - are the moments of C(r). 
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(iv) Gg(0) compensates the negative sum J2x^o Gg (x) and makes Xq positive. Moreover, at 
large £, where £ ~ exp(27r/3), 

G L (0) ~ — - — ~ — (51) 
q[ ' (InO 2 /9 2 ' 1 ' 

giving rise to a positive contact term in the continuum limit. 

In the Table we report some results for G^(0), and the scaling quantities C L (x) (for some 
values of x) and £ 2 Xq ■ 

Among the lattice techniques used to determine Xq m Monte Carlo simulations, there 
are also non-local estimators (see e.g. cooling). Since Xq = "^(Q 2 )) where Q is the total 
topological charge (Q = J2 X Qx) and Q is stable (approximately on the lattice) under local 
changes, in the cooling procedures one changes appropriately the configuration to read Q 
and therefore to determine Xq- But of course this procedure does not leave the two-point 
correlation function invariant: only its zero-mode should be left intact. 

In conclusion we have seen that the continuum limit of the correlation function Gg(x) 
is regular at fixed physical distance and it is given by G q (x). Also its moments Xq,j have a 
regular continuum limit. On the other hand, a singular behavior is found at x = consis- 
tently with reflection positivity and positivity of Xq- These features should also characterize 
the Euclidean correlation function of the topological charge density in the continuum limit 
of lattice QCD. 
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TABLES 



TABLE I. The large- iV limit of NG%(0), C L (x,0) = £ A NG%(x£,0) and £ 2 N X q for 
£ = l,2,4,8,16,oo (corresponding to = 0.41829..., 0.52868..., 0.63901..., 0.74933..., 0.85965..., oo 
respectively). For £ = oo the continuum results are recovered. 



e 


NG%(0) 


C L (l/4,0) 


C L (l/2,0) 


C L (1,0) 


C L (2,0) 


eN X L q 


i 


0.5371 






-0.07888 


-0.001376 


0.199791 


2 


0.2751 




-0.8148 


-0.01634 


-0.001257 


0.171593 


4 


0.1759 


-9.314 


-0.1713 


-0.01615 


-0.001160 


0.162993 


8 


0.1277 


-1.790 


-0.1728 


-0.01547 


-0.001136 


0.160316 


16 


0.09995 


-1.808 


-0.1662 


-0.01535 


-0.001131 


0.159498 


oo 




-1.718 


-0.1651 


-0.01532 


-0.001129 


0.159155 
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FIGURES 
FIG. 1. Plot of In [— C(x)] versus r = \x 
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